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Abstract 

The problem of load balancing in a distribution network under unknown time- 
varying demand and supply is studied. A set of distributed controllers which regulate 
the amount of flow through the edges is designed to guarantee convergence of the solu- 
tion to the steady state solution. The results are then extended to a class of nonlinear 
systems and compared with existing results. Incremental passivity and internal model 
are the main analytical tools. 

1 Introduction 

Cooperative control systems have been widely investigated in a variety of different contexts 
[16, 12, 15, 3]. Less attention has been devoted to cooperative control in the framework of 
dynamical flow networks, with some interesting exceptions [9, 6, 8, 7, 17, 5]. The aim of 
this paper is to study a class of cooperative control algorithms in the context of distribution 
networks under exogenous inputs. 

Main contribution. We analyze and design distributed controllers at the edge which achieve 
load balancing in the presence of time- varying demand and supply (exogenous signals) . The 
role of internal model and incremental passivity is investigated for the problem at hand. 
Similar tools have been used for controlled synchronization and leader-follower formation 
control in e.g. [18, 3, 15, 10] and references therein. We address a different problem and we 
tackle it in a novel way. The load distribution problem is then considered for a more general 
class of systems and this allows us to make a comparison with the results of [2] and [9] . 

Literature review. The literature on the control of flow or distribution networks is wide 
and multi-disciplinary. Here we restrict ourselves to a very small portion of it, focusing on a 
model which takes into account the amount of stored material at the nodes and mass balance. 
This class of systems has been used to model data networks [13] and supply chains [1] for 
instance. Our paper focuses on the problem of stabilizing the flow network to a steady state 
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solution in the presence of exogenous time- varying demand and supply under the scenario in 
which the controllers aim at a uniform distribution of the material among the nodes. This is 
a problem which has attracted considerable attention in the area of parallel and distributed 
computation [16] and has been recently reconsidered for instance in [5] where input and state 
constraints have been taken into account and a connection with [12] has been established. 
The work [5] did not consider the presence of external inputs. A large amount of work on 
the topic of flow control in the presence of disturbances has been carried out in works such 
as [6, 8, 7] where the problem is cast in the robust control framework. The approach in 
our paper is based on the theory of output regulation and to the best of our knowledge this 
has not been considered before. A similar problem has been tackled in [17] but the authors 
restrict themselves to the class of constant disturbances. 

The organization of the paper is as follows. The class of systems under study is introduced 
in Section 2, the design of the edge regulators is carried out in Section 3 and the extension 
to a class of nonhnear system in Section 4. The conclusions are discussed in the last section. 



2 Distribution networks and demand supply balancing 

Consider the system 

X = B\ + Pd (1) 

with X e the state, A G the control vector and d e Mfl, q < n, a disturbance vector. 
The {n x m) matrix B is the incidence matrix of an undirected graph G = {V, E) where 
\V\ — n, \E\ — m. The ends of the edges of G are labeled with a '+' and a Then 

{+1 i is the positive end of k 
— 1 i is the negative end of k 
otherwise 

The system above is a simple model of a flow network [6] and it has been used also to 
model data networks [13] and supply chains [1]. The state Xi G M, i G T := l,2,...,n 
represents the quantity of material stored at the node i, Xk G M, k = 1,2, ... ,m the flow 
through the edge k. The disturbance dj G M represents the inflow or the outflow at some 
node. 

The available measurements are the differences among the quantities stored at the nodes 

namely, z = B^x. 

We assume that each disturbance dj is supposed to be generated by the exosystem 

dj = r^Wj, j = l,...,q, 

where wj G W^^ is the state of the exosystem which describes the evolution of the in- 
flow/outflow j and ^j,Sj are suitable matrices. Considering more general classes of ex- 
osystems is left for future research. We give the system above the compact form 

d = r^'w, 
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where w = (wf . . . wj)^, d = {d^ . . . d^f, S'^ = block.diag(5f, ...,S^),r'^ = block.diag(rf, . . . , T^). 
The model (1) and the overall exosystem (2) return the closed- loop system 

w — S'^w 
X = BX + Pw 
z = B^x 

where by a slight abuse of notation we renamed PV^ simply as P. 

We are interested in the problem of distributing the cumulative imbalance of the network 
due to the in- and out-fiow among the nodes. More formally the problem at hand is as 
follows: 

Load balancing at the nodes Find distributed dynamic feedback control laws 

Vk = ^kVk + J^kZk 

such that, for each initial condition (wq, Xq, r^o), the solution of the closed-loop system (3), 
(17) satisfies lim^^^^ z{t) = 0. 

In what follows we propose a solution to the problem. 



(3) 



(4) 



3 Design of regulators at the edges 

We focus on flow networks whose underlying graph satisfies the following standing assump- 
tion: 

Assumption 1 The graph G is connected. 

The first result concerns the characterization of a "steady state" solution to the problem: 

Lemma 1 Let Assumption 1 hold. For each w solution to w = S'^w, if there exist a function 
Xw : IR+ M"* and a continuously difjerentiahle function x'^ : M+ solution to 

x-^ = BX^ + Pw (5) 

and 

= S^x^ (6) 

then 

llPw 

X ; X^ (7j 

n 

and Xyj — Mw, for some matrix M. If the graph is a tree, then the matrix M is unique. 

Proof: From Assumption 1 and (6), one obtains that for some function 

x'^ : M+ M. Replacing the expression of x^ in (5) one has 

Inx": ^BX^ + Pw, (8) 
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As J\f{B'^) = 7^(l„), with J\f,n the null space and the range of a matrix, multiplying on 
the left both sides of (8) by 1^ yields = as claimed. Replace the latter in (8) to 

obtain YPw = BX^, Y = 1^-^ — In- By Assumption 1 and without loss of generality (up to 
a relabehng of the edges of the graph), BX^ — BaX^ + B^X^ with Ba full-column rank and 
e R"-^ If a solution to YPw = BX^ exists, then A^ = {B^Ba)-'^B^{YPw - B^X^). 
Letting A^ = one obtains A^ = M„w, with M„ = {BlBaY^BlYP. 

If G is a tree, then B is full-column rank and A™ = {B^ B)^^ B^Y Pw . □ 
In what follows, we assume that a solution to (5), (6) exists. Moreover, if m > n — 1, 
then without loss of generality we assume that the first n — 1 columns of B are linearly 
independent and we let the last m — n + 1 components of A*" be identically zero. 

Remcirk 1 Prom (7), by integration, one has 



Observe that x"^ depends on the initial condition and strictly speaking cannot be referred 
to as a steady state solution. Bearing in mind the interpretation of (1) as a flow network 
and of Pw the vector of the inflows and outflows of the network, the integral ^n^(«) ^g 
can be seen as the cumulative imbalance of the network. In other words, if for any given w 
a solution to the load balancing problem exists, then the state at each node equals - up to 
a constant - the cumulative imbalance of the network. 

In the case of a network with no imbalance, i.e. l^Pw{t) = for all t > 0, is a constant 
vector. 

Example 1 Consider the graph depicted in Fig. 1. The graph corresponds to system (1) 
with 

/-I IN 

S = 1-1 , P = 

V 1 -ly 

The solutions of (5)- (6) (with w — d) are as follows 

3 

- A3 H 3 — 
_ ^u; _j_ d\—d2 
2 3 3 ' 

A solution is obtained letting A3 = 0. 

We introduce now a system which generates the control signal A^ in Lemma 1. Consider 
the input A^ associated with the edge /c, with /c = 1, 2, . . . , m. In general, such input may 
depend on all the components of the disturbance vector w. Hence, to generate A^, the 
following system is proposed: 




(9) 

Uk = Hkr]k 



The statement below is immediate. 
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Figure 1: The distribution network considered in the Example 1. 



Lemma 2 For any w solution to w = S'^w, there exists a solution rj^ to (9) such that 
HkVki^) ~ ^ki^) f^^ ^ > 0, where is the kth entry of \^ in Lemma 1. 

Proof: Choose w(0) as the initial condition of (9), then rik{t) = w{t) for all t > 0. As 
= Mw, then it suffice to choose as the kth row of M to have Hkf]'j^{t) = Hkw{t) = 
X^{t) for all t>0. □ 

Remark 2 From the proof of Lemma 1 it turns out that m — n + 1 components of A"' can 
be chosen identically zero. The matrices Hk corresponding to these components are then 
identically zero as well. Hence, for k = n,n + 1, . . . ,m, the system (9) reduces trivially to 
Uk = 0. 

The system (9) is completed by adding control inputs Vki,Vk2 to be designed for guaran- 
teeing that the response of the closed-loop system converges to the desired response for x. 
Hence, we set 

Vk = S'^r]k + Vki /-j^QX 
Uk = Hkr]k + Vk2, A; = 1,2, . . . ,n - 1 

with rjk, Vki G M*^, Vk2 G IR, and Uk = Vk2 for = n, -|- 1, . . . , m. 
We write (10) in the form 

7] = Sr] + vi , . 

X = H7] + V2 ^ ^ 

where rj = {rjj rjj . . . ?7^_i)"^, 5" = ® 5''^, where Cg> denotes the Kronecker product, and 

H=(^^^'y H,= block.diag(i/i, . . . , Hr^.i). 

Observe that by Lemma 2, for any w and provided that vi = 0, V2 = 0, there exists a 
solution rj"^ to (11) which satisfies 



:i2) 



Theorem 1 Consider the system (1), where B is the incidence matrix of a graph G and d 
is a disturbance generated by the system (2). 
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Under Assumption 1, provided that Sj is skew symmetric for each j = 1,2, ... ,q, the 

—T 

dynamic feedback controller (11) with Vi = —H B^x and V2 = —B^x, namely 

fj — Srj — H^B^x /.o\ 
A = Hr,-B^x ^'^^ 

guarantees houndedness of the state of the closed-loop system and asymptotic convergence of 
x{t) to ln(c' + /q ^'^^"'W dg'^ jor some constant d . 

Proof: Consider the overall closed-loop system 

w = S'^w 

X = B(Hr] + V2) + Pw 
77 = Sr] + vi 
z — B^x 

Introduce the new variables x — x — x'^ , r) — r] — rf" . These satisfy 

5 = B{Hri + V2) + Pw - BX" - Pw 
= BHfj + B(Hr)'^ - X"^) + Bv2 
= BHfi + Bv2 

and _ _ 

= Sfi + vi. 

Introduce the Lyapunov function V{x, fj) — ^ {x^x + fj^fj). The function V computed along 
the solutions of system 

X = BHfi + Bv2 /-^^x 
?7 = Sf] + vi 

satisfies V{x, fj) = x'^ {B Hfj-\- Bv2)+fi'^ {Sfj+Vi) . Under the assumption of the skew-simmetry 
of S, one obtains V{x, rj) = x^BHfj + x^Bv2 + fj^vi. Set 

vi = -ifB'^x, V2 = -S^x. (15) 

Observe that by the connectivity of the graph and the definition of x, Vi = —H B^x and 
V2 = —B^x. Then V{x,fi) = —\\B'^x\\'^. Hence, {x,fj) is bounded. By La Salle's invariance 
principle and connectivity of the graph, the solutions to (14) converge to the largest invariant 
set contained in {{x,fj) : B^x = 0} = {{x,fj) : x e 7?.(1„)}. 
Observe that the system (14) with the inputs v as in (15) becomes 

5 = -BB^x + BHfj , , 

r\ = bfj — H B X 
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On this invariant set the system (16) satisfies 

i = BHfi 
fj ^ Sfj 
= B^x. 

Hence, x = Replacing this expression in the equation for x and pre- multiplying both 

sides by 1^, one obtains x^, = 0, that IS IS db constant. Hence x = x — x'^ ^ l„c for some 
constant c. Bearing in mind the expression of x"^ obtained in Lemma 1, then one concludes 
that x{t) l„(c' + Jq ^"•^"'^^^ ds) for some constant c'. □ 

Remark 3 In the case of balanced demand/supply, the state x{t) converges to l„c' for 
some constant c'. Observe that I'^x = 0, that is I'^xit) = l"^a;(0). Hence, l'^x{0) = 

limi^oo = nc' implies that x(t) converges to In^—^^- Hence under the effect of a time- 

varying but balanced demand/supply all the components of the state x{t) asymptotically 
converge to the average of the initial distribution of material at the nodes. 

Bearing in mind the block diagonal nature of the matrices S, H and the definition 
z — B'^x, the dynamic feedback controller (13) can be decomposed as the following set of 
dynamic feedback controllers at the edges: 

Afe = HkTik-Zk, /c = 1,2, . . .,n - 1 

which only requires the knowledge of the difference between the quantities stored at the two 
nodes connected by the edge. As such the proposed controller (13) is fully distributed and 
solves the load balancing problem formulated in Section 2, with $fc = S'^, Ak = —H^, = 
Bk, Tfc = —1. By Remark 2, for /c = n, n -|- 1, . . . , m for which B^ — the edge controller 
becomes a static one, i.e. = —Zk- 

Example 1 (Cont'd) Assume that di — a + /3sm(u!t + tp), with a > /3 > and d2 — a. 
The supply is a periodic fluctuation around a constant value while the demand is a constant. 
Then the matrices S'^ and in (2) write as 

Let Ag = 0. Then, for A; = 1, 2, the matrices Bk which allow to reproduce A^ are 

\ 0), H,= {\ -I 0) 
Then the controllers at the edges 1 and 2 are given by (17) with 5"^ and Bk as above and 

Zi^—Xi^-X2, Z2^—X2+X'i. 

The controller at edge 3 is the static control law A3 = — ^3 — —{xi — xz). 
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Remark 4 (Passivity-based reinterpretation) The proof of Theorem 1 can be reinter- 
preted as follows. In view of Lemma 1, the system 

i ^ B\ 
z — B^x 

is the incremental model associated with system (1). Similarly, by Lemma 2, system 

fj — Sf] + H^v 
u — Hfj, 

where u — u — u'^ and u'^ :— Hrf" , is the incremental model associated with the internal 
model ^ 

f] = St] + H V 
u = Hf] 

The systems are passive with respect to the storage functions Vi(x) = ^x'^x and V2{fi) = \'iffj 
provided that S'^ is skew symmetric. The negative feedback interconnection of the two 
systems, namely 

^ — ^ext ~ 

is passive as well from the input {Xext,Uext) to the output {z,u). The output feedback 

Uext J \ J \ U 

gives asymptotic convergence of the closed-loop system to the largest invariant set where 
z = 0. 

We discuss briefly the difficulties related to the presence of possible state and input 
constraints. 

State constraints. Consider a variation of the model (1) in which the positivity constraint 
on the amount of material stored at the nodes is enforced. The model becomes 

x = {BX + Pw)t 

where (fejA -|- Piw)'^. is the ith component of the vector {BX + Pw)'^ and 



Ci if {xi > 0) or (i = and C,i > 0) 
if {xi = and Q < 0) 



We consider the special case of balanced demand and supply, i.e. l^Pw — 0. As a conse- 
quence, Pw — —BX^ and 

x = i = iBX)t. 
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The function Vi{x) = ^x^x, with x = x — Ix^ and x"^ > 0, satisfies 

V^ix)^x^iB~X)t. 

Observe that 3F{BX)'^ = Yll=i^i{^i^)ti = SF{B\). This shows that the system 

i = {B\)t 

z — B^x 

is passive and the arguments of the previous remark can be used. A formal analysis requires 
to take into account the discontinuity of the system. This is not pursued here for lack of 
space. 

Edge capacity constraints. Constraints on the capacity of the edges can be modeled via a 
saturation function replacing A in (1) with sat(A). Here, sat(A) = (sat(Ai) . . .sat(A^))-^ and 
sat(Afe) = min{|Afc|, c}sign(Afe). Following Lemma 1, let x^,X" be such that 

x"" = Ssat(A^) + Pw 

and M such that sat (A"") = Mw. For the problem to be feasible restrict the set of initial 
conditions wq of the exosystem w — S^w in such a way that ||Mw(t)||oo < c for all t > 
([11]). Consider the incremental model 

i = B[sat(A) - sat (A"")] 
= Ssat(A) - BMw 

To tackle the problem, wc assume the scenario in which at each edge a dynamic observer 
provides w that converges to w asymptotically (or at each edge k there exists an estimator 
which generates a local estimate Wk of w). Consider then the control input 

A = -n{B^x) + Mw 
= -ii{B'^x) + Mw, 

where : R™ — > M'" is a map such that each component is a monotonically increasing 
function which is zero at the origin. The incremental model writes as 

£ = Bsat{-fi{B'^x) + Mw) - BMw 

= fisat(--^(B^x) + Mw + M{w - w)) - BMw. 

The right-hand side is bounded and the solutions exists for all t > 0. Suppose that each 
component of /i is a function whose range is within [— |, |]. Then after a finite time, 
seit{—n{B^x) + Mw + M{w — w)) = —^{B^x) + Mw + M{w — w) and the incremental 
model evolves as 

i = -Bn{B'^x) + BM{w -w). 

Consider the projected state y = Qx, where Q is an (n — 1) x n matrix such that Qln — 0, 
QQ^ = In-i and Q'^Q ^ In - In^n/n. It yields 

y = -QBfi{B^Q^y) + QBM{w-w). 
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The unforced system has a globally asymptotically stable equilibrium^; moreover the forcing 
term is decaying to zero. Since the response of the system is bounded then the state y 
converges to the origin which implies that x converges to 7?.(1„). Then one can proceed as 
in the last part of the proof of Theorem 1. The proposed solution relies on the existence of 
distributed estimators for whose actual design is left as a topic for future research. 



4 Flow networks with nonlinear dynamics at the nodes 

In the previous section, the dynamics describing the evolution of the storage variable at each 
node was given by 

Xi ^ hi\ ^- PiW , i = l,2, (18) 

where hi and pi are the ith row of the incidence matrix B and P respectively. Consider 
now a different case of a flow network in which the way material accumulates at the node is 
described by a non-trivial dynamics, namely 

Xi^ fi{xi) + hi\+PiW, i = l,2, (19) 

with vector of measurements Hi G M™ given by 

Vi = hjxi. 

The nonlinear system (19) allows us to put the results of the paper in a broader context 
and compare them with those in [2], [9] (see the end of the section). Observe that for 
k = 1,2, ... ,m, Uik is either Xi, —Xi or 0. The sum of the outputs over all the nodes 
returns the vector of relative measurements 

n 

z = B'^x = ^ yi. 

i=l 



Each system 

ii = fi{xi) + biX + PiW 
Vi = bJxi, i = l,2, 

is assumed to be incrementally passive. 



(20) 



Assumption 2 There exists a regular^ storage function : R x M x M+ — )■ M+ such that 

dVi dVi X , X, 

'dt d^-^^'^""'^ biX+Piw)+ 
dV ' 

^iUxd + kX' + p,w) < (y, - y'f{X - X'). 



^Take V{y) = then V{y) < and V{y) = is identically zero if and only if B'^Q'^y = 0. This 
implies that y = 0. In fact if this were not true, that is B^Q^y = and y ^ 0, then QBB^Q^y = as well 
and this would contradict that y ^ since QBB^Q^ is a non singular matrix. 

^See [14] for a definition. 
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Remark 5 (A class of incrementally passive systems) Consider the linear dynamics 
at the node (18) and the function Vi = ^{xi — x'j)"^. Then the right-hand side of the inequality 
above becomes 

{xi - x'j){biX + p^w) - {xi - x'j){biX' + piw) 

= {x,-x';)k{\-y) 

= {hJ{x,-x\)f{\-\') 

= {vi - yinx - AO 

which satisfies the dissipation inequality in Assumption 2. 

Suppose that the dynamics fi are equal to VFj, with Fj a twice continuously differentiable 
and concave function. Then the static nonlinearity —fi{xi) is incrementally passive, that is 

(x,-x'Mi^^)-M^'i))<0. 

As a matter of fact fi{xi) — fi{x[) = VFj(xj) — VFi{x'j) = V^Fi{^i){xi — x'^ for some lying 
in the segment connecting Xj, x[. By concavity, V^-Fj(^j) < and therefore {xi — x^){f.j{x,j) — 
fii^'i)) < 0. Hence any system (20) with fi{xi) = VFi{xi) and defined as before satisfies 
Assumption 2. 

Lemma 1 is replaced by the following: 

Lemma 3 For each i = 1,2, . . . ,n, for each w solution to w = S'^w, there exist a function 
: and continuously differentiable bounded functions xf : M+ M that satisfy 



if = Mxf)+b,\''+p,w, i = l,2,...,n 



= E": 



(21) 



i=l 



only if there exists a solution : M+ — )■ M defined for all t > to 



n n 



where f{x) — {fi{x) . . . fn{x)) . If this is the case, then 



A^ 



= 1,2,. 

A- \ _ / MJ{xf) + M2W 



A^ 7 \ 



with A^ G W~, A^ e R"'-"+\ and Mi, M2 suitable matrices. 

Proof: The second equality in (21) and Assumption 1 implies that x'^ = x"^ for all i. 
Replacing the latter in the first equality implies that necessarily x^ must be a solution of 
the inhomogeneous differential equation 

n n 
11 



Suppose that a solution xf exists for all t and let xf = xf for each i. Then the second 
equation in (21) is satisfied by the connectivity of the graph and the properties of the 
incidence matrix. Since xf — x^ for all i, it is seen that the first equation in (21) is satisfied 
if and only if there exists A"' such that 

r[/«) - Pw] = BX\ with r = Mi - /„. 

By connectivity of the graph, the rank of B is n — 1. If the graph has n — 1 edges (i.e. it is a 
tree), then B is full-column rank and, provided that a solution to the previous equation 
exists, it is given by A*" = [B^ B)^^ B^Y [f [x'^) — Pw). If the graph has more than n — 1 
edges, then without loss of generality (up to a relabeling of the edges of the graph) we can 
partition B as {Ba Bi,)^ with Ba full-column rank. Then, provided that a solution to the 
previous equation exists, it is given by A^ = {B'^Ba)~^Bj[Y{f{x'^) — Pw) — BhX]^]. One 
particular solution is obtained for A^ = and A^;' = {B^ Ba)~^BlY{f{x'^) - Pw). □ 

Remcirk 6 If the inflow and outflow are balanced, i.e. = 0, then the solution x'^ to 

(22) exists for all t and is bounded. In fact, consider the system 



y 



n 



and the radially unbounded function V(y) — Then 



n — ' n 

1=1 



By the incremental passivity property of — /«, yfiiy) < for all i and this implies V{y) < 0. 
Hence every solution to the system above is bounded and so is x'^. 

Remark 7 In the case the dynamics at the nodes are all the same, i.e. fi = fj for all i,j, 
then the expression of A"' simplifies as 

This descends from the proof, since by definition of the matrix Y, y/(x^) = 0. 

In the remaining of the section we assume that a solution to (21) exists. 
The parallel interconnection of the n subsystems (20) with input A and output z = ^^Li yi 
returns an incrementally passive systems. Formally 

Lemma 4 The parallel interconnection 

ii = fi{xi) + biX + piw 

in = fn{Xn)+bnX+PnW 



n 
i=l 



h r- 



12 



denoted as 

X = f{x) + BX + Pw 



(23) 



is such that the storage function V{x, x') — X]r=i ^(^j) satisfies 

^if{x) +BX + Pw) + ^{f{x') + BX' + Pw) 

<iz-z')'^iX-X'). 

The proof is straightforward and is omitted. Consider now systems of the form 

Vk = MVk,Vk) ^24) 
Uk = i>k{Vk): A; = l,2,...,n-l, 

with the following two additional properties: 

Assumption 3 For each k = 1,2, . . . ,n — 1, there exists regular functions Wk{r]k, r]'^) such 
that 

Assumption 4 For each k — 1,2, . . . ,n — 1, for each w solution to w — S'^'w, there exists 
a bounded solution rf^ to f]k — (f>k{Vk, 0) such that — ipkiVk)- 

Assume that the system 

■w _ '^nfiVka) I ^n-^Vkb 
Vka — ~ ^ 



n n 

w 
kb 



is forward complete. Initialize the system as ^^^^(O) = a;^(0) and rj^f,{0) = w{0). Then 
Vkait) = ^:it) and 77^,(i) = w{t) for all t > 0. Hence = M^kfiVka) + M^uVkb, k = 
1,2,. . .,n — 1, where M^. and M2k are the kih rows of Mi and M2 respectively. On the 
other hand, = 0, k — n,n -\- 1, . . . ,m. An expression for 0^, A; = 1, 2, . . . , n — 1 is 



n n 
S'^Vkb 



'ipkiVk) = Mikf{r)ka) + M2kr)kb- 

In the special case of nodes with the same dynamics (/j = fj = f for all i,j) ipk{Vk) simplifies 
as M2k'r)kb and a system that satisfies Assumptions 3 and 4 is 

f)k = S% + M^^Vk 
Uk = M2kr]k, 
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with storage function Wk{rik) = ^vlVk- Collect the systems (24) into a system with state 
variable rj = {rji . . . ri^_i)^ , input v = {vi. . . Vm)'^ and output u = {ui . . . UmY i namely 

(25) 



77 = $(77,f) 



with $(77, v) — . . . 0^_i)^, ^(77) = (V'l ■ ■ ■ V'n-i 0^)^. The system is incrementally passive 
from V to u with storage function VF(?7, 77') = X]fe=i Wk{Vk, v'k)- 
The following holds: 

Theorem 2 Let Assumptions 1-4 hold. Suppose that a solution to (21) exists and x'^ is 
hounded. Consider the systems (23), with input A and output z, and (25), with input v and 
output u, interconnected via the relations v — —z + v^xti X — u + Xext- 
The interconnected system is incrementally passive from the input (A^^^ vj^t)'^ to the output 
{z'^ u^Y . Moreover, the feedback {X^xt '^Ixti^ — {—Kz^ 0^)^ , with K a positive definite 
diagonal matrix, guarantees lim(_^+oo z{t) = 0. 

Proof: The feedback interconnection of incrementally passive systems is incrementally 
passive ([14], Lemma 1). Hence 



X = 


f{x) +BX + PW 


z = 




V = 


$(77, v) 


u = 


^(^) 


A = 


U + Xext 


V = 


-Z + Vext 



is incrementally passive from the input (A^^^ '^Ixt)'^ ^'^ output (2;^ -u^)^. The storage 
function U is given by the sum \/ + where V, W are the functions defined above (in 
Lemma 4 and after (25), respectively). 
Let Xext — —Kz, Vext = 0. The system becomes 

X = f{x) + B{^{rj)-KB^x) + Pw 
fj = ^{rj, —z) 
z = B^x 

For a given solution w to w = S'^w, let x^, A"' be as in Lemma 3 and 77"^ as in Assumption 4. 
The functions x^ and 77"" are a solution to the equations above with input (A^^ '^^xtY — ^ 
and output {z'^ u^ f = (O^A""^)^. In fact 

^ f{x'") + B-^{7]'")+Pw 

= fix"") + BX"^ + Pw 

^ $(r^"',o) 

= B'^x,,,. 
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As in [14], by the incremental passivity of the feedback system and the existence of a solution 
[x'^jT]'^) of the feedback system such that z{t) = 0, any other solution {x,r)) with input 
(AL vLf = i-Kz^ 0^)^ satisfies 

V{{x,r]), {x'^,r]'^)) < 

{{z^ u^) - (0^ A-^)) ^ = -z'^Kz. 

Bearing in mind the regularity of U and boundedness of x^, this yields boundcdncss of x. 
In view of the time- varying nature of the system, to infer convergence of z to zero, one can 
resort to Barbalat's lemma. This guarantees convergence under the assumption that z is 
bounded. This in turn requires w bounded, which is the case here since S is skew symmetric. 

□ 

Corollary 1 // (i) fi — f for all i — l,2,...,n, (ii) there exists a twice continuously 
differentiable convex function F{x) such that VF{x) — f{x) and (Hi) l^Pw — for all 
t > 0, then the controllers 

Vk = S'^i^k - M^f^Zk 

Afc = M2k'nk-Zk,, A; = 1,2, . . . ,n - 1, 

and Afe = —Zk, k = n,n + 1, . . . ,m, guarantee limt_>+oo z{t) = 0. 

The closed-loop system given in the corollary above takes the form 

X = VF(x) + BX + Pw, z = B^x 

i) ^ ^T] - M^z, A = M2r] - z , 

where we are assuming that m = n — 1 for the sake of simplicity. This system can be 
compared with similar ones appeared in the recent literature ([2], [9]), where models of the 
form 

X = VF{x) + BX, z = B^x 
ri ^ z, A = -^iv) 

with ip a non-decreasing monotonic non-linearity (such as a saturation function) , were stud- 
ied. The presence of the non-trivial dynamics S in our controller is due to the time varying- 
nature of the external input. In [2], VF{x) has a unique equilibrium at the origin and the 
system x — VF{x) + BX is strictly passive. In [9] it is shown that if the components of 
the vector field VF{x) have different equilibria, VF(x) is strongly concave and ijj introduces 
saturation constraints, then the system's response exhibits state clustering. 

5 Conclusions 

We have presented an internal model approach to the problem of balancing demand and 
supply in a class of distribution networks. Extensions to nonlinear systems have also been 
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discussed. Further research will focus on a detailed investigation of state and input con- 
straints and more complex models of demand and supply. The fulfillment of the internal 
model principle has to be understood for more general classes of nonlinear systems than 
those in Corollary 1. This will shed light on the relation between the results in this paper 
and the saddle-point perspective of [9]. Compared with other papers where the robustness 
to time- varying inputs is studied using a frequency domain approach ([4]), our state space 
approach allows us to consider more general classes of cooperative control systems. 
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